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Compaction of a granular material under cyclic shear 
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Abstract. In this paper we present experimental results concerning the compaction of a granular assembly 
of spheres under periodic shear deformation. The dynamic of the system is slow and continuous when the 
amplitude of the shear is constant, but exhibits rapid evolution of the volume fraction when a sudden 
change in shear amplitude is imposed. This rapid response is shown to be to be uncorrelated with the slow 
compaction process. 

PACS. 45.70.Cc Granular compaction - 61.50.-f Crystalline state 



1 Introduction 



Powder compaction has recently attracted the attention of 
physicists as the prototype of a disordered system without 
thermal fluctuations. For millimetric grains, thermal en- 
ergy is negligible compared to gravitational energy. Once 
poured in a box, the grain packing is trapped in a meta- 
stable configuration and stay in this state unless an ex- 
ternal excitation such as vibrations is imposed. During 
such a process, the particles assembly exhibits a slow and 
complex evolution towards a more compact configuration. 
This evolution has been studied experimentally by Knigth 
et al Q and Nowak et al Q in the case of a vertically vi- 
brated packing of spheres. 

Their experiment consisted in a vertical cylinder full 
of monodispersed beads submitted to successive distinct 
vertical taps of controlled acceleration. The mean volume 
fraction was recorded after each tap. The experiments 
were first performed with taps of constant amplitude jjj. 
They have shown that the increase in volume fraction is 
a very slow process well fitted by the inverse of a loga- 
rithm of the number of taps. The more energetic taps are 
the more efficient for compaction. Nowak et al M have 
studied the compaction under taps of variable amplitude, 
and showed that irreversible processes occur during the 
compaction. Starting from a loose packing, the evolution 
of the volume fraction is not the same when increasing 
the amplitude of vibration as when decreasing. The first 
branch is irreversible, whereas the second is reversible. The 
highest volume fraction is obtained by first increasing the 
tap acceleration then decreasing it back to zero. Recently, 
Josserand et al || have investigated the evolution of the 
packing to sudden change of tap amplitude. They showed 
that the response of the system depends on the history of 
compaction. 



These experimental results have motivated numbers of 
theoretical and numerical studies trying to find a minimal 
model which could exhibit qualitatively similar features. 
Most of the models suggest that the rearrangement of par- 
ticles in a more compact state needs more and more co- 
operative reorganization as the volume fraction increases. 
Absorption-desorption model |^],|],||], geometrically frus- 
trated lattice model random walk model || or ex- 
cluded volume model have been proposed JlO| , |ri| . Most of 
these approaches aspire to describe the vertically vibrated 
system. To our knowledge the question of how the com- 
paction dynamics depends on the method of excitation 
has not yet been raised. In absence of thermal fluctua- 
tions, one can legitimately wonder how the compaction 
dynamic depends on the way the system is excited. 

An alternative way of excitation is horizontal vibra- 
tion. In a previous experiment we have shown that very 
compact and crystalline packings of monodispersed beads 
can be obtained this way |]i"2| whereas vertical vibrations 
only yield random close packings. However, under hori- 
zontal vibration, only the top layers of the packing are af- 
fected and indeed submitted to a periodic horizontal shear. 
The goal of this paper is the study of the response to shear 
of the whole packing. A parallelepiped box full of beads 
is submitted to a horizontal shear through the periodic 
motion of two parallel walls. The resulting compaction is 
analyzed in detail in this paper. 

Deformation of granular material under cyclic solicita- 
tion has been extensively studied in soil mechanics [fl3|| . 
Cyclic loading in a triaxial test shows a compaction but 
no systematic study of the dynamic of compaction for 
long time behavior has been investigated. Cyclic shear 
is also evoked by Scott et al Q in a short note where 
they show that this deformation applied to monodispersed 
beads packing could lead to a partially ordered structure. 
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bottom plate 



Fig. 1. Schematic of the shear cell and position sensor. D is 
the horizontal displacement of the bottom plate and h is the 
height of the packing. Angle of shear is 9 = arctan(Z?//i). 



A shear apparatus has been also developed by Bridgwa- 
ter and Scott Jl5|,[l6| in order to study the segregation 
dynamic during shear. Cyclic shear has been recently ap- 
plied to thermally activated systems such as assembly of 
spherical colloidal particles. A periodic shear deformation 
applied to a colloidal hard spheres glass yields the forma- 
tion of ordered regions |I^,|l8|. However, both the thermal 
excitation and the imposed deformation are in this case 
present. The experiments presented in this paper can be 
seen as the analog of the colloidal experiments at zero 
temperature. The paper is presented as follows. The ex- 
perimental setup is described in section ^[ In section || 
the compaction under constant shear amplitude is inves- 
tigated and the crystalisation process and the slow relax- 
ation are discussed. The section |] reports the study of 
the response of the system to sudden change in the shear 
amplitude. Finally we discuss in section the analogies 
and differences between compaction under cyclic shear 
and vertically vibrated system before giving concluding 
remarks in section ||. 



2 Experimental setup 

The experimental setup is sketched in Fig. [p. The shear 
cell was parallelepiped, and the volume occupied by the 
beads was typically 10.5 cm high, 7.9 cm wide and 10.2 
cm deep. The bottom of the cell was a flat PVC plate at- 
tached to a horizontal linear displacement device driven 
by a stepper motor. The two mobile lateral walls were alu- 
minum plates linked with two hinges to the bottom plate. 
These plates were covered with a thin plastic sheet to in- 
sulate the packing from metallic dust due to particle/wall 
friction. The front and back sides were glass plates fixed 
on the bottom plate. 



The granular packing was confined on the top by a 
rectangular plate mounted on a vertical displacement rail. 
This plate was independent from the rest of the shear cell 
and was free to move vertically only, no horizontal dis- 
placement being allowed. The periodic shear deformation 
was obtained by imposing a periodic horizontal displace- 
ment to the bottom plate, the two mobile lateral walls 
being constrained by two cables to remain in contact with 
the top plate. 

The granular material we used were spherical glass 
beads. Experiments were performed with either monodis- 
persed or bidispersed set of beads. For experiments with 
monodisperse particles, beads were 2.97 ± 0.06 mm in di- 
ameter with a density of 2.52 ± 0.02 g.cm -3 . The cell was 
filled with 1450 g of particles. The bidispersed set of parti- 
cles was constituted of 725 g of the 3 mm beads and 725 g 
of 2 mm beads of the same density. In order to prevent 
the surface deterioration of the beads during long-time 
experiments, beads were coated with silicon oil Rhodorsil 
47V100 (viscosity hundred times viscosity of water). Note 
that the same quantitative results were observed without 
lubrication, but lifetime of the beads was in this case quite 
shorter. 

We measured the volume fraction during the compaction 
process by recording the vertical position of the top plate, 
which went down (resp. up) when compaction (resp. di- 
latation) occurred. Its position was accurately measured 
by a linear position sensor Novotechnik T50 (potentiomet- 
ric transducer). The position resolution was 3 10~ 3 cm cor- 
responding to a 2 10 resolution in volume fraction. Both 
the data acquisition and the horizontal displacement of the 
cell were controlled by the same PC computer. Data could 
be recorded continuously during the shear, or once every 
cycle for the study of the long time evolution. 

The experimental procedure was the following. The 
top plate was removed, the mobile side walls were put the 
vertical position and the particles were poured into the 
cell through a hopper. The mean initial volume fraction 
of the packing obtained by this procedure 4>o = 0.592 ± 
0.008. The top plate was then slowly put into contact with 
the packing. Once the initial random packing was ready, 
the periodic shear deformation was imposed. The lateral 
plates were inclined to an angle +9, followed by an inclina- 
tion to an angle —9. The plates were then put back in the 
vertical position and the volume fraction was recorded. A 
new cycle of shear was then applied and so on. The volume 
fraction could also be continuously recorded. 

During the experiments, we have observed that the in- 
clination angle 9 was the only pertinent control parameter. 
Change in the velocity of the bottom plate (from 1.9 to 
3.8 cm/s) or change in the weight exerted by the top plate 
on the packing (from 1.8 to 3.8 kg) did not affect any of 
the measurements. We also performed experiments using 
asymmetric motion of the lateral plates: the plates were 
oscillating between 9 medmm + 9 and 9 medl u m - 9. In the 
range accessible by our set-up (9 < 12°) no influence of 
the medium position has been observed. The compaction 
process was then only controlled by the amplitude of the 
cyclic shear. 
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Fig. 2. Continuous evolution of volume fraction as a function 
of number of cycles; 9 — 5.4°. 

3 Cyclic shear of constant amplitude 

We have first studied the evolution of a packing submit- 
ted to a periodic shear of constant amplitude 9. A typical 
evolution of the packing volume fraction is presented in 
Fig. ||. In this figure the volume fraction 4> was recorded 
with a resolution of 325 samples/cycle. This plot clearly 
shows a slow and continuous evolution from the initial 
random packing towards higher volume fraction. Super- 
posed to this slow evolution, a quasi-periodic oscillation is 
observed, corresponding to the dilatation occurring dur- 
ing one cycle: the packing is denser in the vertical position 
and alternatively dilates when the walls are inclined at 9 
and —6. This behavior is also observed when the medium 
position is not the vertical. In this case a local maximum 
of the volume fraction is measured at the medium position 
9medium whereas a local minimum is observed when the 
walls are inclined at 9 medmm + 9 and 9 me dium - &■ 

In this paper we do not study the dilatancy properties 
in detail but rather focus on the long time evolution of 
the packing volume fraction. For this purpose, the volume 
fraction is recorded only once a cycle when the walls are 
brought back to the vertical position (stars in Fig. |J). 
A monotonous curve is then obtained, representing the 
evolution of the volume fraction as a function of the 
number of cycles n. The results for three different shear 
amplitudes 8 = 2.7°, 5.4°, 10.7° are presented in Fig. | 
for the case of 3 mm beads. 

The first remark is that the compaction process is a 
very slow process. After 6 10 4 cycles (one week experi- 
ment), the volume fraction still slowly increases. This slow 
compaction appears to be more efficient for large shear 
amplitude than for small shear amplitude as shown by in- 
sert in Fig. |^ representing the volume fraction as a func- 
tion of the shear amplitude 9 after 10 4 (asterisks), 3 10 4 
(crosses), and 6 10 4 cycles (triangles). 

Although curves of Fig. |^ present very small fluctua- 
tions for an individual run, the final volume fraction can 
vary significantly (approximately 10%) from one experi- 



Fig. 3. Compaction under different shear angles : 9 — 2.7° 
(star), 9 = 5.4° (square), 9 — 10.7° (circle). Insert: volume 
fraction of the packing versus shear angle for increasing number 
of cycles: n — 10 4 cycles (asterisks), n — 3 10 4 cycles (crosses), 
n = 6 10 4 cycles (triangles). 



merit to another carried out at the same shear angle. For 
example, Fig. |] shows two compaction curves obtained 
for 9 = 5.4 ° starting from almost the same initial volume 
fraction (0.594 ±0.0008) prepared by the same procedure. 
Evolution of the packing thus seems to be rather sensitive 
to the initial random structure of the grains assembly. In- 
sert of Fig. U presents the same curves in semi-logarithmic 
scale. We have tried to fit our experimental results by the 
inverse logarithmic function proposed by Knight et al for 
the compaction under vertical taps (Eq. 3 of Q). How- 
ever, no reasonable agreement with both the short and 
long time behavior could be found (Fig. |^). Other pro- 
posed fits like the two exponential function pQ] or the 
stretched exponential (jjj and references therein) do not 
provide more convincing agreements. 

A major difference between the periodic shear com- 
paction and the vertically vibrated system is the large 
volume fraction up to 0.693 that can be attained by the 
former. Such a value obtained for monodispersed parti- 
cles is much larger than the volume fraction of a random 
close packing, which is 0.63, according to Scott and Kil- 
gour 2l}| . This means that crystalline arrangements are 
created during the compaction, which perhaps explains 
why the log t behavior is not observed in our system. The 
crystal structure of the packing is indeed observed at the 
wall of the cell as shown in Fig. |5|a. Moreover we checked 
that crystallisation is also present in the bulk by carefully 
removing the particles layer after layer at the end of an 
experiment. The picture in Fig. |^b shows a structure ob- 
served in the bulk. The orientation of the crystal observed 
in the bulk is not always parallel to the walls, suggesting 
that the order is not only wall-induced but nucleates and 
grows in the bulk. 

In contrast, no crystalline alloy structure has been ob- 
served with the bidispersed material made of mixture of 
2 and 3 mm beads. However the evolution of the volume 
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Fig. 4. Compaction curves for 9 — 5.4 ° for two different runs. 
Insert: semi-logarithmic scale. Dotted line is the inverse log- 
arithmic fit <f> = <j)oo — (4>oc — 0o)/[l + Bln(l + n/r)] with 
(^oo = 0.74, 0o = 0.594, B = 0.1 and r = 0.8. 





Fig. 6. Example of angle variation effect during the com- 
paction process. Insert shows a close-up of the first jump (point 
A). 



fraction as a function of number of cycles and shear am- 
plitude presents the same trend using mono or bidispersed 
material. In this case, the maximum volume fraction ob- 
tained for 8 = 10.7° after 6 10 4 cycles is 0.65 to be com- 
pared with 0.693 for the monodispersed packing. 




Fig. 5. Photographs of the 3 mm beads packing, (a) front 
view of the packing, (b) Ordered structure in the bulk of the 
packing (upper layers of the packing were removed). 



4 Response to a change in shear amplitude 

After the study of compaction under constant shear angle, 
we investigated the effects of a sudden change of angle. A 
periodic shear with inclination angle 8\ is first imposed to 
a random packing, and at a given time, the shear ampli- 
tude is suddenly changed to another value 82 — 8\ + A8. 
A typical experiment is presented in Fig. |^. Starting with 
a shear amplitude 8\ — 2.7°, the angle is switched to 
9-2 = 10.7° at n = 5000, and then switched back to 8\ 
at n = 10 4 . As can be seen, increasing the shear angle 
at n = 5000 produces a rapid fall of volume fraction <f> 
followed by a slow and continuous increase. When shear 
angle is decreased back at n = 10 4 , a rapid increase of <fi 
occurs, followed by a slower one. 

Such a response is non-trivial. It first shows that mem- 
ory effects exist in the compaction process. Points A and B 
in Fig. ^| correspond to packings having the same volume 
fraction. However, their responses to the same shear am- 
plitude 8 = 10.7° are different: packing A become looser 
whereas packing B pursues its compaction. This means 
that knowing the volume fraction of the packing is not suf- 
ficient to predict the evolution of the system. Such mem- 
ory effects have been recently put in evidence in vertically 
shaken granular packings 

Another important observation is the rapid character 
of the packing response. Whereas the compaction from 
initial random state is a slow process (Sec. |3|), a sudden 
variation of shear amplitude induces a rapid change of vol- 
ume fraction, which can be regarded as a discontinuity or 
a jump. The amplitude A<f> of this jump can be accurately 
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Fig. 8. Relation between volume fraction jumps and angle 
Fig. 7. Jump amplitude versus age of packing. AO = 5.4° change for aged packings (n > 310 4 ). 
(circles), AO = -5.4° (stars). 



defined and measured with the help of the tangent lines as 
shown on the insert of Fig. ^|. This definition is also valid 
for positive jumps (n = 10 4 on Fig. ^). We have systemat- 
ically studied how the volume fraction jumps depend on 
the different parameters. 

First, we have investigated how the discontinuity A<f> 
is influenced by the age of the packing, namely the time 
elapsed before the angle change is applied. Fig. |?j shows 
the response A<j> to the same angle change Ad = 5.4 ° as a 
function of the elapsed time. The jump \A(/>\ increases with 
age and eventually saturates for n > 2000. The same trend 
is observed for positive jumps (A9 — —5.4°). After 2000 
cycles, the response does not depend on the age anymore. 

Secondly, we have observed that the initial angle does 
not play any role. Changing the shear amplitude from d\ 
to d\ + Ad yield a response Acj) in volume fraction inde- 
pendent of 9\ . In fine, the response of the system for suf- 
ficiently aged packings only depends on the shear change 

Ad. 

Fig. ^ summarizes jump measurements carried out for 
the number of cycles n > 3 10 4 . Data points align along a 
straight line crossing the origin: 

A(j) = -aAd (1) 

with a = 1.45 10~ 3 and AO in degrees. In conclusion, the 
response of the packing to sudden change amplitude ap- 
pears to be simpler than first thought. The rapid variation 
of the volume fraction induced by the change is simply 
proportional and opposite to the angle change. 

5 Discussion 

Two kinds of evolution are observed in the system of cyclic 
sheared packing: a slow and continuous compaction when 
the shear amplitude is constant and a rapid response when 
a change in shear amplitude is imposed. The striking fea- 
ture is that the rapid dynamic is uncorrelated with the 



slow compaction: the variation of the volume fraction A<f> 
induced by a sudden change A9 in shear amplitude is inde- 
pendent of the state (j> of the system, and depends linearly 
on Ad. This observation suggests to split the evolution of 
the volume fraction in two terms: 

4>{n) = <t> s iow(ri) - ad. (2) 

The first term 4> a i ow (n) corresponds to a slow and continu- 
ous evolution and the second term is proportional to angle 
d. This description is compatible with Fig. |8[ when a sud- 
den change in 9 is imposed, one observes a corresponding 
jump ctAd in volume fraction, <j) s iow being a continuous 
function. 

One way of interpreting these results is to say that the 
periodic shear deformation plays two roles. On one hand, 
it allows for the particles to slowly and continuously rear- 
range to form a more compact packing. It thus profoundly 
affects the structure of the packing in an irreversible way, 
by introducing more and more order. This contribution is 
4>siow On the other hand, it introduces an additional dis- 
order superimposed to the ordered structure, represented 
in Eq. || by the negative contribution —ad. These two con- 
tributions are uncorrelated. In this framework, the slow 
and rapid dynamics would then correspond to changes in 
the packing structure that are of complete different na- 
ture. 

However, such interpretation remains speculative as 
long as we are not able to relate the macroscopic behav- 
ior of the volume fraction to local rearrangement of the 
particles. An important question is the discrimination be- 
tween the rearrangements involved in the slow compaction 
process from those involved during the rapid compaction 
or dilatation subsequent to a sudden change in shear am- 
plitude. Some numerical simulations suggest the existence 
of collective motion of cluster and of individual motion of 
particles which could perhaps be related to the slow and 
rapid dynamics 

At this stage, it is interesting to compare the results 
presented here about the compaction of a granular pack- 
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Fig. 9. Compaction under variable shear angle. 8 is increased 
linearly from 2.7° to 10.7° then decreased to 2.7° and finally 
increased again up to 10.7°. Rate of increase or decrease is 
±2.7 10 -4 "/cycle. The dashed line has a slope equals to —a. 



6 Conclusion 

In this paper we have presented experiments about the 
compaction of a granular media under cyclic shear. First 
we have put in evidence that this method of compaction 
leads to crystals when the material is made of monodis- 
persed particles. By studying the response of the system 
to sudden change in shear amplitude we have shown that 
discontinuities in volume fraction are observed, with am- 
plitudes proportional to the angle change. The evolution 
of the packing is thus composed of two dynamics: a slow, 
continuous and monotonous compaction and a rapid evo- 
lution occurring during the shear jump. This observation 
suggests that the rearrangement of the particles in the 
packing associated to the two dynamics are of different 
nature. However, the next important step would be to ex- 
perimentally relate the compaction process with the evo- 
lution of the internal structure of the packing. 

This work has benefited from fruitful discussions with R. Monas- 
son, A. Lemaitre and S. Luding. It would not have been pos- 
sible without the technical assistance of F. Ratouchniak. 



ing under cyclic shear with the results obtained previously 
in the case of a vertically vibrated packing. In order to do 
so, we have carried out the annealing experiment analog 
to the one described in Q for the tap excitation: the am- 
plitude of shear is continuously increased then decreased 
and increased again. The evolution of the volume fraction 
is plotted in Fig. ^|. As observed in the tap experiment, 
one get an irreversible branch and a reversible one. 

However, in our case, the reversible branch is simply 
a straight line parallel to -ad. In the light of the above 
interpretation about the response of the system to sud- 
den shear changes, this behavior can be explained as fol- 
lows: slowly increasing the shear amplitude yields a quasi- 
saturated state, where 4> s i ow (n) does no longer vary with 
n. By then changing the angle 8, one explores the linear 
variation of </> with 8. The slight shift between the de- 
creasing and increasing branches in the reversible regime 
comes from the slight evolution of <j> s iow(n). Hence this 
shift increases (resp. decreases) when decreasing (resp. in- 
creasing) the rate of change in angle. 

There thus exist several similarities between the two 
modes of compaction. However, some differences have to 
be noted. First, cyclic shear induces the crystallisation of 
an assembly of monodispersed beads, whereas such order 
is not reported in the experiments of vertical taps. The 
small aspect ratio (10 particles in a cross-section) of the 
tap experiment could perhaps inhibit the crystalisation. 
A second difference is the amplitude of the fluctuations. 
According to [|| vertical taps induce fluctuations from one 
cycle to another of the order of 0.8% of the mean volume 
fraction, whereas in our shear experiment the variation of 
volume fraction is less than 0.1%. Periodic shear is per- 
haps a less violent way of exploring the different packing 
configurations. 
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